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Abstract  -  An  integrated  thermal,  micromagnetic,  spin-momentum-transfer  (SMT)  model 
was  developed  to  study  the  effect  of  SMT  on  the  programming  current  required  for 
thermally  assisted  magnetic  random  access  memory  (MRAM).  The  thermal  portion  of 
the  model  is  used  to  compute  Joule  heating  by  the  spin-polarized  current,  and  it  is  based 
on  a  Crank-Nicolson  inhomogeneous  heat  equation  solver.  The  magnetic  portion  of  the 
model  is  based  on  a  micromagnetic  Langevin  dynamic  Landau-Lifshitz-Gilbert  solver 
including  SMT  torque.  Simulations  of  thermally  assisted  magnetization  reversal  of  0.09 
pin  MRAM  elements,  heated  by  passing  current  through  the  barrier  separating  the  pinned 
and  free  layers  were  performed.  The  free  layer  of  the  MRAM  elements  was  switched 
using  a  magnetic  field  at  fixed  heating-SMT  current  bias.  Results  suggest  that  a  spin- 
polarized  heating  current  can  be  used  to  lower  the  programming  current  required  to  write 
thermally  assisted  MRAM  if  the  direction  of  the  heating  current  is  properly  synchronized 
with  the  reversal  field. 
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One  difficulty  in  developing  high-density  magnetic  random  access  memory 


(MRAM)  lies  in  developing  a  memory  element  with  long  data  retention  time  and  low 
error  rates  that  can  be  programmed  with  currents  that  will  not  damage  the  MRAM 
element  and  can  be  generated  in  a  practical  semiconductor  device.  These  problems 
result  in  part  from  an  increasing  sensitivity  to  thennal  activation  as  MRAM  is  scaled  to 
smaller  dimensions  reducing  the  energy  required  to  reverse  the  memory  element.  In 
order  to  improve  stability  as  MRAM  element  dimensions  are  reduced,  increases  in 
magnetic  anisotropy  or  magnetic  thickness  are  required.  Unfortunately,  at  small  device 
dimensions  these  methods  increase  programming  fields  and  currents  to  impractical  levels. 

Methods  for  writing  MRAM  using  heat  in  combination  with  applied  magnetic 
fields  or  spin  momentum  transfer  (SMT)  have  previously  been  described  to  address  this 
issue.  [1]  Promising  heat-field  or  magnetothennal  (MT)  MRAM  methods  involve  the 
development  of  new  materials  with  high  spin-polarization  and  low  Curie  temperature  or 
the  use  of  an  additional  low  Neel  temperature  antiferromagnet  to  pin  the  free  layer. 
SMT-heat  methods  are  also  promising,  but  work  needs  to  be  done  to  develop  stable 
devices  that  can  be  written  at  safe  current  densities.  In  spite  of  the  fact  that  the  heat-SMT 
and  MT-MRAM  devices  have  many  similarities  in  their  structures,  little  or  no  work  has 
been  performed  to  combine  the  techniques  to  further  reduce  programming  current.  This 
may  in  part  be  because  the  SMT  critical  current,  Ic,  for  reversing  the  magnetization  of  a 
mono-domain  thin  film  ferromagnet  by  SMT  in  zero  temperature  simulations  is  only 
weakly  dependent  on  magnetic  field  (H)  until  H  is  very  near  the  uniaxial  anisotropy  field 
of  the  film,  H|<.  [2]  In  spite  of  this,  it  is  worth  studying  the  combined  effect  of  H  and 
SMT  at  elevated  temperature,  since  Joule  heating  in  MT-MRAM  devices  requires  current 
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8  2* 

densities  close  to  those  required  for  SMT  reversal  (~10  A/cm“),  thin  ferromagnetic  films 
at  -100  mn  dimensions  are  not  strictly  mono-domain,  and  thennal  fluctuations  in  real 
magnetic  layers  cause  reversal  at  H<Hk  on  short  time  scales.  In  order  to  study  the 
possibility  of  an  efficient  combined  SMT-MT-MRAM,  a  heat  equation  solver  was 
developed  and  integrated  with  a  Langevin  dynamic  (LD)  micromagnetic  solver. 
Simulations  of  a  130  x90  x  3.6  mn  elliptical  nickel-iron  MRAM  cell  are  presented  herein 
to  demonstrate  the  model. 

Micromagnetic  models  for  studying  SMT  at  non-zero  temperature  exist  are  often 
based  on  the  stochastic  Landau-Lifshitz-Gilbert  (sLLG)  equation  combined  with  an 
additional  SMT  torque  term: 

^  =  -ym  x  Heff  -  y am  x  (m  x  // ^ )  -  y J  SI  m  x  (m  x  m  p ) .  (1) 

Here,  the  first  and  second  terms  on  the  right  hand  side  are  the  usual  precession  and 
damping  terms,  and  the  third  term  represents  SMT  from  a  fixed  ferromagnetic  layer 
oriented  in  direction  mp.  “f  ’  is  the  time,  y  the  gyromagnetic  constant,  a  the  damping 
parameter,  m  a  unit  vector  representing  the  direction  of  magnetization, 
jsi  =  (r\hi)/(2eM,V),  T|  the  spin  polarization,  Ms  the  saturation  magnetization,  V  the 

computational  cell  volume,  and  Heff  =  Hexc+Hd+Hk+Hext+Hi+Hns,  the  total  local  vector 
magnetic  field  within  the  ferromagnetic  object  due  to  exchange,  demagnetization, 
anisotropy,  external  sources,  current,  and  thennal  agitation. 

Many  of  the  details  of  the  LD  micromagnetic  solver  algorithm  used  in  this  work 
have  been  previously  described.  [3]  The  Gaussian  distributed  Hns  field  sequence  is 
calculated  using  the  Box-Muller  algorithm,  with  a  standard  deviation  of 
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Hns  =  .j2akBT /(yMsVAt)  and  zero  mean.  Hi  is  computed  for  all  current  flowing  within  the 
mesh  using 


BJRles<) 
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In  this  notation,  “test”  indicates  the  point  where  the  B  field  from  a  current  filament  of 
length  L  =  ^ L2X  +  L2y  +  L2_  at  the  “source”  point  is  computed.  Thus, 


Kst  =  Xtes,X  +  ytes,y  +  ZteJ  aIld  ^source  =  XsourceX  +  V  source 9  +  2  source2  •  The  deltaS  are 


Ax:  =  x,  ,  —  x  ,  Av  =  v.  ,  —  v  ,  and  Az  =  z.  ,  —  z 

test  source  ’  J  J  test  s  source  ’  test  source 

The  mx  (in  x  m  )  representation  of  the  SMT  torque  is  probably  the  most 

commonly  accepted  form.  [4]  Addition  of  this  tenn  to  equation  1  has  an  interesting 
consequence  in  that  it  produces  a  current  and  magnetization  orientation  dependent  change 
in  the  effective  damping  parameter  of  the  material.  [5]  Consider  the  case  where  mp  is 
parallel  to  Heff.  In  this  case,  the  sLLG  equation  can  be  rewritten  as 

=  -yrh  x  Heff  -  yam  x  (m  x  Heff )  with  tt  =  a  ±  J SI  j H eff  .  SMT  torque  changes  the 

effective  damping  causing  magnetic  excitations  to  damp  out  more  slowly  or  quickly  in 
time,  depending  on  the  orientation  of  mp.  The  a  value  in  the  thermal  Hns  term  is  not 
rescaled.  The  difference  between  the  damping  in  the  deterministic  and  stochastic 
portions  of  the  SMT  modified  sLLG  equation  implies  that  SMT  torque  should  change  the 

effective  magnetic  temperature  of  the  sample  to  T  =  (a!  a)T ,  when  the  spin  polarization 
of  the  current  is  collinear  with  the  magnetization  of  the  free  layer.  Using  switching 
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measurements,  it  is  difficult  to  separate  T  from  a  lowering  of  the  energy  barrier,  AE, 
since  the  measurements  are  dependent  on  the  ratio  AE/keT.  In  any  case,  the  effective 
increase  in  T  or  lowering  of  AE,  should  facilitate  magnetic  reversal  at  lower  H  or  T. 

This  work  thus  focuses  on  reducing  the  programming  current  of  MT-MRAM  and  through 
the  use  of  SMT. 

Heat  flow  due  to  Joule  heating  at  the  tunnel  barrier  is  simulated  using  the  heat 
equation  modified  for  use  in  nonunifonn  materials  and  to  include  temperature  dependent 
material  parameters: 

Pm  (x>  y  ,z)Cp  (x,  y,z)dr^y,Z)  =  v  ■  [k(x,  y,z  )VT(x,  y,z)]+  g(x,  y,z,  T,  t)  (3) 

Here,  pm  is  the  local  density,  Cp  the  local  heat  capacity,  T  the  local  temperature,  k  the 
local  thermal  conductivity,  and  g(x,y,z,T,t)  is  the  time  and  temperature  dependent  local 
power  generation  in  the  simulated  device.  For  Joule  heating  in  materials  with 
temperature  dependent  resistivity,  the  power  generation  term  is  expressed  as 

g(x,  y,  z,  T,  t)  =  J2(x,  y,  z,  t)p(x,  y,  z,  T ) ,  (4) 

where  the  local  current  density,  J,  is  dependent  on  time,  and  the  local  resistivity,  p,  is 
dependent  on  temperature.  Equation  3  was  discretized  using  a  Crank-Nicolson  finite- 
difference  approach.  The  resulting  set  of  matrix  equations  are  efficiently  solved  in  3D 
using  a  sparse-matrix  preconditioned  successive-over-relaxation  solver 

The  integrated  thermal-micromagnetic-SMT  simulation  method  is  outlined  in 
Figure  1.  Note  that  the  block  labeled  “RG  Scaling”  indicates  that  renormalization  group 
scaling  of  the  micromagnetic  simulation  input  parameters  is  applied  to  account  for  errors 
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resulting  from  truncation  of  short  wavelength  excitations  in  the  simulated  spin  lattice  due 
to  the  coarse  grained  mesh  (5  x  5  x  3.6  nm3).  [6] 

A  cross-section  of  the  simulated  MRAM  element  is  shown  in  Figure  2.  The 
device  sandwiches  a  free  layer  between  two  oppositely  oriented  pinned  layers.  The  free 
layer  is  separated  from  one  pinned  layer  by  an  AlOx  (tunnel  barrier)  spacer  and  from  the 
other  by  a  Cu  spacer.  Since  the  AlOx  barrier  magnetoresistance  is  much  greater  than  that 
of  the  Cu  spacer,  this  structure  can  have  a  large  magnetoresistance.  The  oppositely 
oriented  pinned  layers  effectively  double  the  spin  polarization  in  the  free  layer  since 
across  one  spacer,  spin  polarized  electrons  are  reflected  and  those  across  the  other  are 
transmitted.  In  order  to  simulate  this  stack,  an  effective  r|=0.9  was  thus  used.  Thermal 
parameters  and  geometry  used  for  simulation  are  listed  in  Table  1.  The  free  layer  is 
composed  of  permalloy,  Ms  =  8xl05  A/m,  A  =  le-1 1  J/m,  a  =  0.02,  and  Tc  =  850  K.  The 
tunnel  barrier  has  a  resistivity  of  ~5  T>|xm2,  which  assuming  a  200  mV  reliability 
threshold,  safely  permits  about  300  pA  of  current. 

Joule  heating  results  are  shown  in  Figure  3.  The  free  layer  temperature  is  unifonn 
to  better  than  1%.  The  average  temperature  is  thus  used  for  the  micromagnetic 
simulation.  Micromagnetic  simulation  results  are  shown  in  Figure  4.  The  combination 
of  applied  field,  SMT,  and  temperature  sufficiently  lower  AE  so  that  the  sense  layer 
magnetization  reverses  within  the  write  cycle.  Figure  4(A)  shows  that  increased 
temperature  strengthens  the  dependence  of  the  SMT  Ic  on  H  and  decreases  the  SMT  Ic 
value.  Figures  4(B)-4(D)  show  that  at  low  Iwrt  neither  H  nor  H  combined  with  SMT  can 
reverse  the  bit.  At  intermediate  Iwrt,  only  SMT  -  H  can  program  the  bit,  and  at  higher  Iwrt, 
H  alone  can  be  used  to  program  the  bit.  For  the  present  work,  tunneling  hotspots  through 
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the  barrier  were  also  considered.  Simulations  restricting  tunnel  current  to  an  area  of 


about  10%  of  the  AlOx  spacer  show  more  than  a  50%  reduction  in  SMT  Ic. 

LD  micromagnetic  simulation  suggests  that  increasing  temperature  increases  the 
dependence  of  SMT  on  applied  field.  A  method  for  writing  an  MRAM  element  utilizing 
a  combination  of  SMT,  heat,  and  applied  field  has  thus  been  suggested.  SMT  and  applied 
magnetic  field  allow  thermally  assisted  writing  of  an  MRAM  element  without  the  need  to 
heat  the  free  layer  above  its  Tc  at  currents  lower  than  needed  if  only  Joule  heating  and  H 
were  used  to  overcome  AE.  Switching  speed  was  found  to  range  from  5  ns  at  327  pA  to 
greater  than  50  ns  at  150  pA.  For  the  particular  geometry  and  set  of  materials  studied, 
the  use  of  SMT  provides  roughly  a  15%  reduction  in  write  current.  Simulation  including 
tunneling  hotspots  at  10%  of  the  total  junction  area  shows  switching  at  less  than  150  pA, 
and  can  potentially  reduce  the  write  current  by  50%  below  the  uniform  current  simulation 
results.  Further  simulations  have  shown  that  decreasing  the  Tc  of  the  free  layer  to  500  K 
allows  a  device  of  the  same  dimensions  to  be  built  that  could  be  written  at  less  than  100 
pA  through  the  memory  element,  in  spite  of  the  fact  that  100  pA  only  provides  about  50 
K  of  heating. 
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Figures 


Figure  1.  Integrated  micromagnetic-thermal-SMT  simulation  method. 

Figure  2.  Simulated  MRAM  element.  TV  =  thermal  via,  Ih  =  current  used  to  produce 
the  write  field,  FM-clad  =  ferromagnetic  cladding,  Iwrt  =  write  current,  and  Ir(j  =  read 
current.  The  bit  is  written  by  passing  current  Ih  through  the  write  conductor  and  tapping 
some  current  through  the  bottom  TV  to  produce  SMT  and  heat.  The  direction  of  Iwrt  must 
be  synchronized  with  Ih. 

Figure  3.  (a)  Temperature  as  a  function  of  time  and  current,  and  the  input  write  current 
waveform,  (b)  Temperature  distribution  in  the  plane  of  the  free  layer  20  ns  after  applying 
175  pA.  The  contours  are  spaced  10  K  apart,  and  the  average  temperature  of  the  free 
layer  is  456.5  K±2K. 

Figure  4.  (A)  Average  SMT  switching  current  computed  at  fixed  H  and  T  with  T|  =  0.45. 
(B)-(D)  show  the  benefit  of  including  SMT  plus  Joule  heating  in  the  MT-MRAM  device 
at  write  currents  ranging  from  150  pA  (B)  to  200  pA  (D).  SMT  reduces  the  write  current 
by  about  25  pA.  Inset  of  (B)  shows  SMT  switching  at  150  pA  with  10%  area  tunneling 
hotspot. 

Table  1.  Thermal  simulation  parameters  for  the  model  MRAM  bit.  Simulated  thermal 
resistance  between  FM-FREE  and  CMOS  is  3.6e6  K/W.  The  substrate  and  top  conductor 
are  assumed  to  be  a  constant  heat  bath  at  300.  Neumann  boundaries  are  assigned  to  the 
vertical  edges  of  the  simulated  volume. 
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